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ABSTRACT 

We exploit XMM-Newton archival data in a study of the extended X-ray emission 
emanating from the Galactic Centre (GC) region. XMM-Newton EPIC-pn and EPIC- 
MOS observations, with a total exposure time approaching 0.5 and 1 Ms respectively, 
were used to create mosaiced images of a 100 pc x 100 pc region centred on Sgr 
A* in four bands covering the 2-10 keV energy range. We have also constructed a 
set of narrow-band images corresponding to the neutral iron fluorescence line (Fe i 
Ka) at 6.4 keV and the K-shell lines at 6.7 keV and 6.9 keV from helium-like (Fe 
XXV Ka) and hydrogenic (Fe xxvi Lya) iron ions. We use a combination of spatial 
and spectral information to decompose the GC X-ray emission into three distinct 
components. These comprise: firstly the emission from hard X-ray emitting unresolved 
point sources; secondly the reflected continuum and fluorescent line emission from 
dense molecular material; and, thirdly, the soft diffuse emission from thermal plasma 
in the temperature range, kT k, 0.8-1.5 keV. 

We show that the unresolved-source component accounts for the bulk of the 6.7- 
keV and 6.9-keV line emission and also makes a major contribution to both the 6.4-keV 
line emission and the 7.2-10 keV continuum flux. We fit the observed X-ray surface 
brightness distribution with an empirical 2-d model, which we then compare with a 
prediction based on an NIR-derived 3-d mass model for the old stellar population in 
the GC. The X-ray surface brightness falls-off more rapidly with angular offset from 
Sgr A* than the mass-model prediction. One interpretation is that the 2-10 keV X-ray 
emissivity increases from « 5 x 10^'' erg s~^ M©"^ at 20' up to almost twice this value 
at 2'. Alternatively, some refinement of the mass model may be required, although it 
is unclear whether this applies to the Nuclear Stellar Cluster, the Nuclear Stellar Disc, 
or a combination of both components. 

The unresolved hard X-ray emitting source population, on the basis of spectral 
comparisons, is most likely dominated by magnetic cataclysmic variables, primarily 
intermediate polars. We use the X-ray observations to set constraints on the number 
density of such sources at angular offsets between 2' — 20' from Sgr A* (projected 
distances at the GC of 4.6-46 pc). Our analysis does not support the conjecture that 
a significant fraction of the hard X-ray emission from the GC originates in very-hot 
(~ 7.5 keV) diffuse thermal plasma. 
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1 INTRODUCTION 

The centre of our Galaxy represents an intriguing labo- 
ratory for the study of high energy processes and inter- 
actions. Within the central 100 pc, the rise in the stel- 
lar mass density is complemented by a concentration of 
gaseous, largely molecular matter, in a region known as 
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the Central Molecular Zone jMorris fc Serabvn|[T99^ '). Re- 
cent star formation coupled with supernova explosions ev- 
ery few thousand years serve to inject high-energy parti- 
cles into the mix, which power no n-thermal emission and 
drive outflows l|Crocker et al.|[201ll : and references therein). 
At its heart, our Galaxy hosts a currently-quiescent super- 
massive black hole, Sgr A*, with a mass of 4 x lO'' M0 
jSchodel et al.ll200^ : iGhez et al.ll2008l '). Quite plausibly re- 
cent activity on Sgr A* may also have energised the sur- 
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roundins resion dSu nyaev fc Churazov 199^: 



199e ; iMurakami et al. 2000; IPonti et al.ll2010l : ICapelli et al 



Kovama et al 



2013) • Given this context, it is not surprising that, even 



when the contribution of highly-luminous X-ray binaries 
and transients is excluded, the central region of our Galaxy 
exhibits a rich variety of bo th compact X-ray s ources and 



diffuse X-ray emission ( e.g., Muno et al 



2004al: |Park et all 



l2004l : IWang et al.ll2006l : iMuno et al.ll2008l . |2009| ) . 

One focus of research over the last forty years has 
been the nature of the unresolved hard X-ray emission, 
known as the Galactic Ridge X-ray Emission (GRXE), 
which can be traced along the Galacti c Plane, reaching a 
maxi m um surface brig h tness at the GC (IWorrall fc Marshalll 



19831: 



1 19851 : iKovama et al.l Il986l . 1 19891 : 



.Warwick et al ^ 

Revnivtsev et al.ll2006l ). This unresolved X-ray 



emission is 



known to exhibit many K-shell emission lines from highly- 
ionised atoms, with the most prominent of these lines be- 
ing He-hke (Fe xxv Kq) and H-hke (Fe xxvi Lya) iron 
at 6.7 and 6.9 keV respectively. The flux ratio of Fe xxv 
Ka/Fe xxvi Lya suggests the origin is in optically-thin 
thermal p lasma emissi on at a tem perature ^ 7.5 keV (e.g., 
Kovama et al. 1986. : Kovama et al. .2007a|: lYamauchi et al.l 



20091 ). If this emission is diffuse in nature, its temperature is 



too high for it to remain gravitationally bound to the Galac- 
tic Plane and it w ould instead escape with a v elocity exceed- 
ing 1 000 km s"^ (|Worrall Marshalll ITgS^ : iKovama et all 
1 19861 ). The kinetic energy required to sustain the plasma 
(~ 10*^ erg s~^), is equivalent to the occurre nce of at least 
one supernova explosion ev ery thirty years (|Kaneda et al.l 
Il997l : lYamasaki et al.lll997l '). This far exceeds current esti- 
mates of the supernova rate within the region. Furthermore, 
the spectra of supernova remnants (SNRs) exhibit thermal 
plasma components with temperatures typically in the range 
of kT ~ 0.2-0.8 keV and seldom, if ever, with kT > 3 keV. 
These difficulties call into question whether the emission can 
be truly diffuse in nature. 

Early Chandra deep-field observations revealed a large 
popul ation of faint point sources with in the GC re- 
gion (Wang et al."2002l: iMuno et all 12003'). More recently, 
iRevnivtsev et al. (200it ) have reported that within a partic- 
ular deep-field observation (/n, &ii = +0.08° , —1.42°), more 
than 80 per cent of the continuum and line emission 
around 6.7 keV can be resolved into point sources. A high 
concentration of discrete sources is perhaps to be expected, 
given that the central 2° x 0.8° region o f the Milky Way 
conta ins ~ 1 p er cent of the Galactic mass (|Launhardt et af] 
120021 '). In fact, iMuno et all (|2006l ') compared the spatial dis- 
tribution of the X-ray emission from point sources with the 
stellar mass distr i bution model derived from the NIR map of 
iLaunhardt et al.l (|2002l ') and, within the limits of the statis- 
tics, found a good correlation between them. They conclude 
that the majority of point sources with Lx ^5 lO^** erg s~^ 
are likely to be magnetic cataclysmic variables (CVs). In 
such systems, an accretion shock is generated above the 
white dwarf surface, which heats the accreting material to 
temperatures kT > 15 keV. The resulting highly-ionised, 
optically-thin plasma cools and, eventually, settles onto the 
white dwarf. The resulting X-ray spectrum comprises a com- 
plex blend of contributions generated at differing tempera- 
tures, densities an d optical depths fe.g.. lCropper et al1ll999l : 
lYuasa et al.ll2O10l ). As well as a hard continuum, the X-ray 
spectra of magnetic CVs are characterised by a complex of 



iron-K lines including strong He-like, H-like and fiuorescent 
components, the latter ar ising from relatively co ld iron in 
states from Fe I to xvil l|Hellier fc Mukaill2004l ). Support 
for the conjecture that magnetic CVs represent an impor- 
tant GC source population also comes from studies of the 
hard X-ray emitting, low-lumino sity X-ray source po pula- 
tion found locally in our Galaxy ijSazonov et al.ll2006l ): and 
the distribution of the GRXE across the cent ral quadrants 
of t he Galactic Plane (IRevnivtsev et al.ll200^ . 

lUchiv ama et ahi (|201ll ) recently investigated the Fe xxv 
Ka emission from both the Galactic Ridge and GC regions 
using Suzaku data and compared the spatial profiles of the 
line emissi on with the stellar mass distribution models con- 
sidered bv lMuno et al. 

liOM)- 

They concluded that the dis- 
tribution of the line emission along the Galactic Ridge is 
consistent with an unresolved point source origin but ar- 
gued that an additional component is required to account 
for the excess emission seen within the GC region. In con- 
trast, based on a detailed spectral ana lysis of Suzaku ob - 
servations in the Galactic bulge region, lYuasa et al.l (|2012l ) 
concluded that virtually all of the unresolved X-ray emis- 
sion could be explained in terms of the integrated emission 
of two point source components, magnetic accreting white 
dwarfs, which give rise to the hard emission extending from 
5-50 keV, and stellar coronal sources which might account 
for the softer {kT 1.2-1.5 keV) thermal emission evident in 
the bulge fields. 

In this paper, we report on a detailed investigation of 
the spatial and spectral characteristics of the X-ray emission 
observed within the central 100 parsecs of the Galaxy based 
on data from XMM-Newton. The emphasis of our study is 
on constraining the model for the unresolved, hard X-ray 
emitting sources and setting limits on what fraction of the 
emission might reasonably be described as diffuse in nature. 
In the next section, we give details of the XMM-Newton 
observations selected for this analysis and the procedures 
followed in the data reduction. We also describe the tech- 
niques employed to construct mosaiced X-ray images cover- 
ing a range of broad and narrow energy bands and to ex- 
tract representative X-ray spectra. In §3 we outline how the 
spatial and spectral information can be modelled in terms 
of three distinct emission components. We then discuss in 
§4 the contribution that unresolved low-luminosity sources 
make to the GC X-ray emission and whether, in light of our 
results, there is any requirement for truly diffuse, very-hot 
thermal plasma emission in the GC. Finally §5 gives a brief 
summary of our conclusions. Throughout thi s work, the dis- 
tance to the GC is assumed to be 8 kpc ijGillessen et al.l 
120091 ). 



2 OBSERVATIONS AND DATA REDUCTION 

The XMM-Newton Science Archive (XSA) was used to se- 
lect observations, publicly available as of 2011 July, with 
pointing positions within 0.7° of (^n, &ll) = (0°, 0°). A mini- 
mum effective expos ure of 1000 s in at least one of the three 
EPIC instruments (jStriider et al.ll200ll : pTurner et al.ll200ll ) 
was set as a requirement. Similarly only observations with at 
least one of the EPIC cameras operating in either full-frame 
(EE) or extended full- frame (EFF) mode were considered. 
The majority of the observations were carried out using the 
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medium filter, altliougli eitlier tlie tliin or thick filter were 
deployed in some instances. Table [T] provides a summary of 
the set of XMM-Newton observations used in the current 
analysis and related information. 

The data reduction was based on SAS vll.O (hereafter, 
SAS). For every observation listed in Table O the Observa- 
tion Data Files for each of the EPIC instruments, pn, MOS-1 
and MOS-2, were reprocessed using the SAS tasks epchain 
and EMPROC. Each dataset was screened for periods of soft 
proton flaring through the creation of full-field light curves 
for the energy range 10-12 keV (utilising only single-pixel 
events). Good Time Interval (GTI) selections were made 
by setting a suitable threshold for each instrument, chosen 
to filter out those periods when the instrument background 
count rate was highly variable and clearly in excess of the 
local quiescent level. The resulting exposure times for each 
instrument after the application of the GTI filtering are re- 
ported in Table [1] The total "live" time summed over the 
set of 50 pn observations and the 56 MOS-l/MOS-2 obser- 
vations was, respectively, ~ 0.5 and ~ 0.95 Ms. The event 
lists were filtered so as to include only single- and double- 
pixel events for pn (PATTERN<=4) and up to quadruple-pixel 
events for MOS-1 and MOS-2 (PATTERN<=12). Out-of-Time 
(OoT) events have a non-negligible impact on pn data and 
appropriate corrections were applied. As a preliminary step, 
the SAS task epchain was used to generate OoT event lists, 
filtered in an identical way to the primary datasets. 

2.1 Construction of the image mosaics 

Images from each observation and from each camera were 
constructed in four "broad" bands covering the 2-10 keV 
bandpass (see Table[2|. Instrument exposure maps were also 
made for each of these bands. In the case of the pn data, 
events with energies in the range 7.8-8.3 keV (which en- 
compasses the prominent copper Ka instrumental line) were 
excluded when constructing the very-hard band image. Sim- 
ilarly for the MOS very-hard band, events in the 9.0-10.0 
keV range were discarded (since the MOS response rapidly 
declines in this energy range). 

Images were also made in a number of narrow bands 
encompassing prominent Fe-K spectral lines. Table [2] lists 
the energy ranges used in the present study corresponding 
to the neutral (or near-neutral) iron fluorescence line (Fe I 
Ka) at 6.4 keV and the K-shell lines at 6.7 keV and 6.9 keV 
from He-like (Fe xxv Ka) and hydrogenic (Fe xxvi Lya) 
ions of iron. Two narrow bands sampling the "continuum" at 
energies immediately below and above the iron-line complex 
were also utilised (see Fig. [l}. 

The first correction applied to the pn images involved 
the construction of a set of OoT images (one per band) from 
the OoT event lists. These OoT images were scaled to the 
estimated OoT event fraction (6.3 and 2.3 per cent in FF 
and EFF modes respectively) and subtracted from the raw 
pn image using the ftools task farith (ftools v6.11, 
iBlackburnlfiggg ). 

The interactions of high-energy cosmic-ray particles, 
with the EPIC detectors are known to give rise to an in- 
strumental background consis ting of both a hard continuum 
and fluorescent line features (jLumb et al.ll20()^ ). A correc- 
tion for the instrumental background was applied to both 
the pn and MOS images by making use of EPIC filter wheel 



Table 1. Details of the XMM-Newton observations used in this 
work. All observations were made with a Medium filter in place 
unless otherwise stated. 



ObsID Observation date Exposure time (s) 

MOS-1 MOS-2 pn 



Total time: 942464 910071 487036 



t Thick filter deployed, 
t Thin filter deployed. 



0112970701 


2000- 


-09- 


-11 


23110 


23116 


16951 


0112970401 


2000- 


-09- 


-19 


22871 


22671 


13712 


0112970501 


2000- 


-09- 


-21 


14650 


14758 


4716 


0112971001 


2000- 


-09- 


-24 


8500 


8898 


7356* 


0112971501 


2001- 


-04- 


-01 


6205 


7242 


2929 


0112971901 


2001- 


-04- 


-01 


5215 


5804 


1339 


0112972101 


2001- 


-09- 


-04 


22500 


22515 


16447 


0111350101 


2002- 


-02- 


-26 


38944 


39164 


31139* 


0111350301 


2002- 


-10- 


-03 


8074 


8097 


5914* 


0202670501 


2004- 


-03- 


-28 


26554 


29254 


9064 


0202670601 


2004- 


-03- 


-30 


34440 


37048 


24267 


0202670701 


2004- 


-08- 


-31 


81716 


84000 


38114 


0202670801 


2004- 


-09- 


-02 


102082 


106098 


62475 


0302882501 


2006- 


-02- 


-27 


8167 


8407 


4595 


0302882601 


2006- 


-02- 


-27 


2858 


3153 


966 


0302882701 


2006- 


-02- 


-27 


4892 


5097 




0302882801 


2006- 


-02- 


-27 


7167 


7072 


2231 


0302882901 


2006- 


-02- 


-27 


6667 


6575 


1876 


0302883001 


2006- 


-02- 


-27 


6764 


6268 




0302883101 


2006- 


-02- 


-27 


10491 


10101 


1593 


0302883201 


2006- 


-03- 


-29 


6132 


6138 


2541 


0302883901 


2006- 


-09- 


-08 


6380 


6285 


2366 


0302884001 


2006- 


-09- 


-08 


6377 


6389 


3319 


0302884101 


2006- 


-09- 


-08 


6074 


5683 


1225 


0302884201 


2006- 


-09- 


-08 


6480 


6086 


4363 


0302884301 


2006- 


-09- 


-09 


6379 


6181 


4014 


0302884401 


2006- 


-09- 


-09 


5440 


5145 


2748 


0302884501 


2006- 


-09- 


-09 


8154 


7763 


5065 


0506291201 


2007- 


-02- 


-27 


24438 


19010 




0402430701 


2007- 


-03- 


-30 


30390 


29680 


8827 


0402430301 


2007- 


-04- 


-01 


52528 


42381 


16728 


0402430401 


2007- 


-04- 


-03 


36809 


31887 


15841 


0504940101 


2007- 


-09- 


-06 


6323 


5918 


2831 


0504940201 


2007- 


-09- 


-06 


9696 


9588 




0504940401 


2007 


-09- 


-06 


6117 


6227 


4396 


0504940501 


2007 


-09- 


-06 


5912 


5929 


4050 


0504940601 


2007 


-09- 


-06 


3552 


3061 




0504940701 


2007- 


-09- 


-06 


6419 


5726 


3355 


0511000101 


2008- 


-03- 


-03 


1582* 


1781* 




0511000301 


2008- 


-03- 


-03 


5431* 


5238* 


606* 


0511000501 


2008- 


-03- 


-04 


6020* 


6229* 


1528* 


0511000701 


2008- 


-03- 


-04 


6422* 


5828* 


2825* 


0511000901 


2008- 


-03- 


-04 


6028* 


4828* 


4303* 


0511001101 


2008- 


-03- 


-04 


6518* 


6529* 


4352* 


0511001301 


2008- 


-03- 


-04 


4730* 


2950* 


2342* 


0505870301 


2008- 


-03- 


-09 


12275 


12073 


4704 


0505670101 


2008- 


-03- 


-23 


79223 


73233 


44097 


0511000201 


2008- 


-09- 


-23 


6319* 


6325* 


3856* 


0511000401 


2008- 


-09- 


-23 


3688* 


2481* 


3682* 


0511000601 


2008- 


-09- 


-23 


6518* 


6328* 


4198* 


0511000801 


2008- 


-09- 


-27 


6409* 


6224* 


1931* 


0511001001 


2008- 


-09- 


-27 


6518* 


6529* 


4352* 


0511001201 


2008- 


-09- 


-27 


6518* 


3522* 


722* 


0554750401 


2009- 


-04- 


-01 


28142 


20968 


24703 


0554750501 


2009- 


-04- 


-03 


41733 


39752 


27885 


0554750601 


2009- 


-04- 


-05 


36523 


34838 


23601 
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Table 2. Energy ranges (in keV) of the four broad bands covering 
the full 2 — 10 keV bandpass and of the narrow bands encompass- 
ing the three spectral lines of iron. 



Table 3. Factors (/) used in interpolating between the 
continuum-low and continuum-high narrow band measurements 
to give the continuum underlying the line image. 



Broad 


bands 




Medium 


2.0 - 


4.5 


Hard 


4.5 - 


6.0 


Iron 


6.0 - 


7.2 


Very-hard 


7.2 - 


10.0^ 


Narrow 


bands 




Continuum-low 


5.70 - 


- 6.15 


Fe I Ka 


6.270 


- 6.510 


Fe XXV Ka 


6.525 


- 6.825 


Fe XXVI Lya 


6.840 


- 7.110 


Continuum- high 


7.20 - 


- 7.65 



^The 7.8 — 8.3 keV spectral region encompassing the Cu Ko in- 
strumental line was excluded in the case of the pn data, whereas 
for the MOS data the very-hard band was truncated at 9.0 keV. 





Line 


pn 


MOS 


Fe 


6.4-keV 


0.72 


0.64 


Fe 


6.7-keV 


0.53 


0.435 


Fe 


6.9-keV 


0.325 


0.245 



Table 4. Percentage spillover from one iron-line band to another 
as a result of the limited spectral resolution of the EPIC cameras. 







pn 


MOS 


Fe 6.4-keV to Fe 


S.7-keV 


8.5 


2.7 


Fe 6.4-keV to Fe 


S.9-keV 


5.0 


4.1 


Fe 6.7-keV to Fe 


S.4-keV 


2.2 


1.9 



o 





1 1 
Low 




1 

6.4 


1 

6.7 


1 

6.9 




High 


























6 7 
Energy (keV) 



Figure 1. The full-field spectrum from the pn camera (ObsID: 
0202670801) in the region of the Fe-K lines. The underlying con- 
tinuum was sampled in two narrow bands (marked "Low" and 
"High") at energies immediately below and above the Fe-K com- 
plex. The various narrow bands are delineated by the vertical 
lines. 

closed (FWC) data obtained from the background analysis 
website maintained by ESAQ. With the sky signal blocked, 
the FWC exposures provide a relatively clean measure of 
the instrumental background and, therefore, may be used 
to model and subtract the instrument background. Merged 
event files containing all available FWC data were selected 
for both the pn FF (Rev. 266-Rev. 2027) and pn EFF (Rev. 
355-Rev. 1905) modes. Merged FWC event files for the MOS 
cameras in FF mode were also obtained with the caveat that 
the FWC data for MOS-1 were split into two categories: be- 
fore (Rev. 231-Rev. 961) and after (Rev. 961-Rev. 2027) the 
failure of CCD-6 on 2005 March 9. The MOS-2 FWC event 
list contained data taken from Rev. 230-Rev. 2027. Using 
the SKYCAST software, the FWC data (initially in detec- 

^ http : //xinm2 . esac. esa. int/ external/ximn_sw_cal/background/ 
f ilter_closed/index. shtml. 



tor coordinates) were cast into images in the sky-coordinate 
frame applicable to the particular observation. The number 
of counts measured in the corners of the detector (i.e., the 
regions of the detector shielded from the sky) in the actual 
observation were then divided by the equivalent number of 
counts in the FWC image to give a FWC scale factor for each 
band, instrument, and observation. A correction for the in- 
strument background was then applied to each observation 
by subtracting the scaled FWC image. 

As a final step prior to the combining of the images into 
a mosaic, the images were inspected visually for undesirable 
artefacts such as scattered fight features and bright tran- 
sient sources. These artefacts were spatially masked when 
encountered and a corresponding correction applied to the 
exposure maps. 

In the production of the image mosaics, observation 
0202670801 (centred at RA, Dec of 266.4413°, -29.028972°) 
was used as a reference field. In the case of the pn cam- 
era, the image data deriving from all the other observations 
were resampled and co-added to the reference image. This 
process was repeated for each band and also for each expo- 
sure map. Flat-fielding then involved dividing each co-added 
image by the appropriate co-added exposure map. For the 
narrow-band images we used the co-added 6.0-7.2 keV ex- 
posure map in this fiat-fielding step. The same process was 
used for the MOS datasets, although in this case the MOS-1 
and MOS-2 image data and exposure maps were combined 
prior to the flat-fielding step. The resulting mosaiced images 
have a format of 648 x 648 pixels (4" pixels) corresponding 
to a 43.2' X 43.2' field. At the distance of the GC this field 
size covers a projected extent of 100 x 100 pc. 

The upper-left panel of Fig.gjshows the 2.0-10.0 keV pn 
mosaic produced by summing the four broad-band images. 
The corresponding exposure map is shown in the upper-right 
panel. Fig. [5] (lower-left panel) also presents a blow-up of the 
region to the north-east of Sgr A*, where there is extended 
high surface brightness X-ray emission. 

In order to trace the distribution of the iron-line emis- 
sion, it is necessary to make a correction to the narrow-band 
images for the in-band continuum emission underlying the 
emission lines. Here the approach was to calculate an inter- 
polated "continuum" image by multiplying the narrow-band 
continuum- low and -high images (see Table (2)1 by a factor / 
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Figure 2. Image mosaics based on XMM-Newton observations covering a central 43.2' X 43.2' field (100 X 100 pc at the GC). The 
images are constructed on a celestial coordinate grid with east to the left. Upper-left panel: 2—10 keV pn image. The adjacent colour 
bar is calibrated in units of pn counts per pixel per 20 ks. Upper-right panel: The corresponding pn exposure map. The adjacent colour 
bar is calibrated in units of logio of the exposure (sec). Lower-left panel: A blow-up of the 2-10 keV image showing the region to the 
north-east of Sgr A*. Lower-right panel: The 2—10 keV pn image overlaid by a spatial mask which removes the most luminous sources 
(see ^2.2| I. The longitudinal and latitudinal spatial cuts are represented by the two intersecting white rectangles, which have long-axis 
dimensions of ±22' and full widths of 8' (see i|2.3l l. The red boxes (of dimension 8' X 9') define the four regions from which spectral data 
were extracted (see i|2.4l l. In three of the panels the white dotted lines represent a Galactic coordinate grid with a 15' spacing. 



and 1 — / respectively, and then summing the two products. 
Full-field spectra extracted from a number of the GC obser- 
vations (see Fig. [TJ were used to estimate appropriate values 
of / for the three Fe-lines, for both the pn and MOS channels 
(see Table |3]). Continuum-subtracted line images were pro- 
duced by subtracting the resultant interpolated continuum 
images from the line images (with the former scaled so as 
to account for the different bandpasses of the narrow-band 
line and continuum images. 

The need for a further correction to the line images 
arises due to the limited spectral resolution of the pn and 
MOS cameras. This results in some spillover of signal from a 
given iron line into the adjacent line band or bands (see Fig. 
[T]). Here we used the X-ray spectral-fitting package (xSPEC, 
V12.6.0) to estimate the spillover factors for spectra domi- 
nated either by refiection or hot-thermal emission. The for- 



mer may be modelled as a bright Fe I Ka line at 6.4 keV 
plus an associated Fe i K/3 line at 7.11 keV (with 11 per cent 
of the Kof normalisation) superimposed on a power-law con- 
tinuum with spectral slope F ~ 1.7. Our simulations showed 
that the two reflection lines contaminate the 6.7-keV and 
6.9-keV narrow-bands as summarised in the first two lines 
of Table |4] (separate values are quoted for the pn and MOS 
instruments since their spectral resolution and energy re- 
sponse differ). Similarly, the spillover of the 6.7-keV He-like 
Fe Ka emission into the 6.4-keV narrow-band was simulated 
using was a solar-abundance apec model with the pl asma 
temperature set at kT = 7.5 keV jSmith et al.ll200ll '): the 
results are given in the third line in Table ID Final Fe-line 
images were produced by subtracting appropriate fractions 
of the other line images. Hereafter, we will often refer to the 
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Fe I Ka, Fe xxv Ka and Fe xxvi Lya lines simply as the 
6.4-keV, 6.7-keV and 6.9-keV lines, respectively. 

Fig. [3] shows the pn continuum-subtracted images for 
the three iron lines. The lower-right panel of Fig. [3] also 
shows the very-hard (7.2-10.0 keV) band image. Note that 
all of these images have been heavily smoothed with a cir- 
cular Gaussian function of width cr = 5 pixels. The bright- 
est 6.4-keV line emission is concentrated to the north-east 
of S gr A* in a region with a high molecular cloud den- 
sity (iKovama et al.il996l : iPark et al.ll2004 iPonti et al.ll2010l : 
ICapelli at al.ll2011 ). In contrast, the surface brightness in the 
6.7-keV and 6.9-keV iron-line images show a smooth increase 
towards Sgr A*, as might be expected if these lines arise pre- 
dominan tly in the integrated emi s sion of unresolved point 
sources l|Muno et all |2006| . l2009l : iRevnivtsev et al.l l2007l : 
iHong et al.ll2009D . On the other hand, the very-hard band 
image can be readily interpreted in terms of a combination 
of the two emission distributions noted previously. 



2.2 Spatial masking of point sources 

The presence of luminous point sources within the field of 
view introduces a degree of confusion noise into the im- 
ages, which can be removed via spatial masking. In order 
to produce such a spatial mask, a catalogue of sources was 
constructed based on the Second XMM-New ton Serendipi- 
tous Source Catalogue (the 2XMMi catalogue: I Watson et al.l 
[2009. ). Specifically, 2XMMi sources within the GC field with 
2.0-12.0 keV flux greater than ~ 1.5 x 10"^^ erg s"^ cm"^ 
were included in a preliminary source list. The mosaiced 
images were then inspected at the source locations and a 
number of possibly spurious sources discarded. Conversely, 
a small number of sources not included in the 2XMMi cata- 
logue, but clearly visible in the images, were added manually 
to the list. The result was a catalogue of 84 sources. A spatial 
mask was then made by excising circular regions of radius 
25" at each source position (corresponding to the removal 
of roughly 80 per cent of the source response). Larger ex- 
cision circles were used for two very bright sources to the 
north-east of Sgr A*, for the Arches cluster and for the pul- 
sar wind nebula, GO. 13-0. 11. We also excluded the brightest 
regions of the Sgr A East complex (using a contour based on 
its 6.7-keV Fe-line emission). The full mask array was then 
applied to every image employed in our study prior to any 
smoothing; the lower-right panel of Fig. [2] shows the result 
in the case of the 2-10 keV pn image. 



2.3 Spatial cuts 

Our investigation of the various components of the diffuse 
X-ray emission present in the GC requires both a spatial and 
spectral modelling approach. We based our spatial modelling 
on cuts through the various images in Galactic longitude 
and latitude centred on the position of Sgr A*. The white 
rectangles in Fig. [2] (1 ower-right panel) indicate the regions 
of the mosaiced images which were thereby sampled. The 
cuts extend out to offset angles of ±22' with a ±4' extent 
perpendicular to the cut axis. The spatial resolution along 
each cut is 1'. We excluded segments of the cuts for which 
more than half the area was already excluded by the source 
mask. This resulted in the loss of the central ±2' region 



(which is confused by Sgr A East) in both cuts and the 
eastern-most 4' section of the longitudinal cut. The results 
from this analysis are presented in 



2.4 Spectral extraction 

X-ray spectra were extracted from four regions within the 
central 100 pc zone, corresponding to the sky areas sampled 
by the spatial cuts at offsets from Sgr A* of between 4' and 
13' (see the lower-right panel of Fig. [5]). The dimensions 
of each region are therefore 8' x 9', with the longer axis 
orientated along the cut direction. Luminous point sources 
were excluded from the extracted regions corresponding to 
those masked out in ^2.2\ 

For simplicity, we extracted the spectra from a single 
pn observation (0202670801) for which the exposure time 
after GTI-filtering was ~ 62 ks. The wide extent of the dif- 
fuse emission in the GC field makes the selection of a back- 
ground region for use in spectral analysis problematic. We 
opted, therefore, to determine the background spectra from 
the appropriate FWC dataset. FWC spectra were selected 
from geometrical areas on the detector matching those of 
the four sky regions. A simple scaling of the FWC spectrum 
so as to match the sky exposure then proved adequate as 
the instrumental background. Subsequent spectral analysis 
then requires all the relevant GC emission components as 
well as the cosmic X-ray background (CXB) to be included 
in the spectral model (see H3.4p . The SAS tools ARFGEN and 
RMFGEN were also used to produce the appropriate auxil- 
iary response and response matrix files. The spectra were 
grouped so that each contained at least ten counts per bin. 



3 SPATIAL AND SPECTRAL MODELLING OF 
THE GC EMISSION 

3.1 Distinct emission components 

A number of unresolved X-ray emitting components are evi- 
dent in the GC region, which can be distinguished by virtue 
of their distinctive spatial and spectral characteristics. Here, 
for simplicity, we attempt to model this emission in terms 
of just three major components. 

The first is the emission associated with the unresolved 
hard X-ray emitting point sources. In the present analysis 
we have masked out the sources brighter than 1.5 x 10~^^ erg 

cm~'^ in the 2-12 keV band which, at the distance of the 
GC, corresponds to a 2-10 keV X-ray luminosity threshold 
of ~ 10"^"^ erg s~^. Hence here we are dealing with the inte- 
grated emission of sources less luminous than this. Our spec- 
tral description of this component is in terms of a hard con- 
tinuum (taken to be thermal bremsstrahlung with kT = 7.5 
keV) with associated Fe-line emission at 6.4, 6.7 and 6.9 keV 
(i.e. neutral. He-like and H-like iron Ka emission). We as- 
sume that unresolved sources give rise to the bright central 
concentration of emission, which is most clearly discerned in 
the 6.7-keV and 6.9-keV images (Fig. [3]). 

The second major emission component observed in the 
GC is that associated with dense molecular clouds. The most 
obvious signature is the iro n fluorescence line at 6.4 keV 
iKovama et al.l l2007bl . l2009l : iPark et all l2004l : iPonti et all 
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Figure 3. The distribution of iron-line emission and the associated very-hard continuum in the central 100 pc of the Milky Way. Upper- 
left panel: the neutral iron fluorescence line at 6.4 keV. Upper-right panel: the He-like iron line (Fe XXV Ka) at 6.7 keV. Lower-left panel: 
the H-like iron line (Fe xxvi Lya) at 6.9 keV. Lower-right panel: the very-hard (7.2—10 keV) band image. All four images have been 
smoothed with a Gaussian mask function of width a = 5 pixels. The adjacent colour bars are calibrated in units of pn counts per pixel 
per 20 ks. The white dotted lines represent a Galactic coordinate grid with a 15' spacing. 



l2010l : [Capelli et al.ll2012l ). This fluorescence is very likely ex- 
cited by hard X-ray illumination of the clouds. The Thomson 
scattering of the incident X-ray flux back into our line-of- 
sight gives rise to a second spectral characteristic, namely 
the presence of a hard power-law X-ray continuum (with 
photon index F in the range 1.5-2.0). Dense molecular ma- 
terial is known to be very unevenly distributed in the GC 
with a concentrat ion to the north-east of Sgr A* {e.g., 
iTsub oi et al.l Il999l ). This explains the highly asymmetric 
surface brightness distribution seen in the 6.4-keV image, 
which is also reproduced in the very-hard band image (see 
Fig.©. 

A third emitting component is that arising from diffuse 
thermal plasma, most likely energised by recent supernova 
explosions and, in some localised regions, by the colliding 
stellar winds of massive stars. This relatively soft emission 
component is pervasive throughout the GC region; however, 
the highest surface brightness region outside of the central 



few arcmins, where emission associated wi th the Central 
Cluster and the Sgr A Ea st SNR dominates (|Baganoff et al.l 
I2OO3I : iMuno et al.il2004bl ). is again located to the north-east 
of Sgr A*, where the density of molecular material is high. 
The temperature and, possibly, metal abundance of this soft 
thermal plasma may vary with location. In modelling this 
component we assume that the thermal structure can be ad- 
equately represented by two temperature components, one 
at kT « 0.8 keV and the other at kT x 1.5 keV. The former 
accounts for much of the line emission evident in the X-ray 
spectrum of the GC region below ~ 4 keV, the most promi- 
nent spectral line being the Ka line of He-like sulphur at 
2.4 keV. The 1.5-keV component, whilst also contributing 
significantly to the spectrum below 4 keV, is hot enough to 
produce 6.7-keV iron-line emission. 

This paper focuses on an investigation of the spatial 
distribution and spectral characteristics of the emission as- 
sociated with the unresolved point source population. An 
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Figure 4. Spatial cuts through the GC images from the pn camera. The left-hand panels show the longitudinal cuts (where -|-ve offset 
is to the Galactic east of Sgr A*) and the right-hand panels show the latitudinal cuts (where +ve offset is to the Galactic north of Sgr 
A*). The spatial resolution is 1' along the cut direction. The cuts have a full width of 8' in the perpendicular direction. The black points 
and histograms represent the data. The best-fitting unresolved source model is shown as the red histogram. The data segments used to 
fit the source component are identified by the vertical dotted lines. Upper panels: Cuts through the 6.7-keV iron-line image. The blue 
histogram represents the residual 6.7-keV emission once the source component is subtracted (only -|-ve residuals are considered); this 
excess signal may be associated with the kT ^ 1.5 keV thermal emission component. Upper-mid panels: Cuts through the 6.9-keV iron- 
line image. Lower-mid panels: Cuts through the 6.4-keV iron-line image. The orange histogram represents the residual 6.4-keV emission 
once the contribution of the unresolved-sources is subtracted (only -|-ve residuals are considered); this excess signal is associated with 
the fluorescence of dense molecular clouds. Lower panels: Cuts through the 7.2-10 keV band image. The orange histogram represents 
the residual emission once the contribution of the unresolved-sources is subtracted (only +ve residuals are considered); this excess signal 
is associated with reflection from dense molecular clouds. 
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analysis of the softer thermal plasma emission will follow in 
a later publication (Heard & Warwick, in preparation). 

3.2 Spatial distribution of the unresolved sources 

As a means of investigating the spatial distribution of the 
unresolved source population, we have employed an empir- 
ical 2-d surface brightness model which can be matched to 
the available spatial information. We adopt a model in which 
the surface brightness, S, decreases as a power-law function 
of the angular offset 9 (measured in arcmin) from Sgr A*. 
The observations demonstrate very clearly that the rate of 
decrease is faster perpendicular to the Galactic Plane than 
along the plane. We accommodate this by including an ex- 
ponential fall-off with the latitudinal angle where (j) is 
measured with respect to the Galactic latitude of Sgr A*. 
The resulting model was: 



S = Ne'''exp{-^) + Sfg, 



(1) 



where the parameters A'^, a, and (f)sc, represent the normal- 
isation of the power-law component, the power-law index, 
and latitudinal scale height (in arcmin), respectively. The 
surface brightness is presumed to peak at the position of 
Sgr A*, although in the model we actually constrained the 
power-law component interior to S = 1' to the value at 
1'. The model also includes a constant foreground surface 
brightness Sfg. We set the latter equal to 10 per cent of the 
value of the variable component aX 6 = 12', cji = Q' . The 10 
per cent factor matches the decrease in surface brightness 
in the 6.7-keV Fe xxv Ka line measured by Suzaku as the 
offset from Sgr A* increases from 12' to 1° I Uchivama et al.l 
I2OIII ). In the present analysis Sfg represents the contribu- 
tion of sources outside of the central few hundred parsec, 
i.e., sources located ei t her in the Galactic bu lge or Galactic 
disc l|Hong et al.ll2009l : lUchivama et al.ll2"oTlh . 

3.3 Modelling the spatial cuts 

Longitudinal and latitudinal spatial cuts derived from the 
6.7-keV, 6.9-keV, and 6.4-keV iron-line images and the very- 
hard 7.2-10 keV continuum image are shown in Fig. U (in 
all cases these are measurements from the pn camera) . Here- 
after, we use the shorthand: Long -|-ve; Long -ve; Lat -|-ve; 
and Lat -ve to identify the different segments of these cuts, 
as per the labelling in Fig. 2] 

A rapid decline in the surface brightness as the distance 
from Sgr A* increases is evident in most cases (with the 
notable exception of the 6.4 keV and very-hard continuum 
measurements in the Long -|-ve region), with the rate of 
decline being somewhat faster in the latitudinal direction. In 
order to match our empirical 2-d surface brightness model 
to these data, we derived a set of model distributions in the 
format of the 1-d cuts, corresponding to an a versus (psc 
parameter grid. Here the process was: (i) construct a model 
image in our standard format for the specified a and ^sc; 
(ii) apply the source spatial mask, thereby mimicking the 
filtering of the actual data; and (iii) determine the spatial 
cut distribution (as per ^2.3[) . 

We then selected the 6.7-keV and 6.9-keV line data per- 
taining to the two "cleanest" regions, namely the Long -ve 
and Lat -|-ve segments in the offset range 2'-16', for fitting. 




Scale height 0,^ (arcmin) 

Figure 5. The variation of as a function of the source model 
parameters, a and t/iac. The contours represent the 68, 90, and 
99 per cent confidence intervals. The results from the pn data are 
shown in red and those from the combined-MOS data in blue. 
The dashed vertical line is drawn at 18.2'; this corresponds to 
the scale height of 45 p c reported for the Nuclear Stellar Disc by 
iLaunhardt et alj 1I2OO2I') (see text). 



The fitting process involved the determination of four free 
parameters, a, 4'sc and two normalisations (one for each 
line) via ^^-minimisation. The results for the two "interest- 
ing" parameters, a and 4>ac are shown in the form of a 2-d 
confidence contour plot in Fig. [S] The best-fit values were 
q=0.601q q3 and ^ac^lS.Gl]^ ^ arcmin, respectively. The 
for the best-fit was 109.0 for 54 degrees of freedom (dof), re- 
fiecting the fact that there is some (modest) scatter over and 
above the predicted level of the noise. Nevertheless, the best- 
fit model does account for the bulk of the 6.7-keV and 6.9- 
keV line emission in the Long -ve and Lat -|-ve regions (see 
Fig.|4]). As a check we have repeated exactly the same anal- 
ysis for the combined-MOS data. The derived constraints 
on the two model parameters a and 4'sc from the MOS are 
also shown in Fig. [5] the MOS resuhs (a=0.56lo;o2 and 
(;/>sc=17.0lJ;4 arcmin; = 79.6 for 54 dof) are fully con- 
sistent with those obtained from the pn channel, indicating 
that our modelling process is robust. 

As a further test, using the best-fit values for a and 
(fisc derived from the pn data, have also we fitted the Long 
-ve and Lat -|-ve cuts (2'-16') for the 6.4-keV iron line and 
very-hard band data, which entailed determining the best-fit 
model normalisation for the respective bandfl The results, 
shown in Fig. UJ indicate that that the unresolved source 
model also accounts for the bulk of the emission measured 
in these bands within the fitted regions. 

Finally the models derived for the Long -ve and Lat -|-ve 
regions were also extrapolated to the Long -|-ve and Lat -ve 
regions (with all the parameters including the normalisa- 
tions fixed). Fig. U shows the results. In the Long -|-ve seg- 
ment, the unresolved source model generally under-predicts 



^ In the case of the 7.2-10 keV band, we obtained our best result 
when we subtracted a constant offset from the cut data of 7.1 X 
10~^ photon/cm^/s/arcmin^ - this may reflect the limitations of 
the background subtraction in this continuum band. 



10 V. Heard and R. S. Warwick 



the observed signal. This is perhaps not surprising given 
that this is the region to the north-east of Sgr A*, where X- 
ray fluorescence and soft thermal emission are very promi- 
nent (Fig. [2] & [Sjl . However, the model does provide quite 
a reasonable description of the 6.9-keV line emission in this 
region. 

In the Lat -ve segment, there appears to be an enhance- 
ment in 6.7-keV line emission over and above the source 
model contribution. This excess extends out to an offset of 
~ fO' (~ 25 pc). There is some evidence for excess 6.4- 
keV emission in this same region, but in this case extending 
only to ~ 5'. Interestingly, Chandra data provide evidence 
of some unusual emission clumps to the south of Sgr A*, 
which arguably form part of a larger-sca le bipolar feature 
l|Baganoff et al.l l2003l : fMorris et al.ll2003h . Plausibly, these 
are all elements of an extended structure linking back to the 
Sgr A East region (Heard & Warwick, 2012 , in prepara- 
tion). Nevertheless, the source model again provides a good 
match to the 6.9-keV and very-hard band cuts in the Lat 
-ve region. 

In summary, our modelling of the GC spatial cuts leads 
us to the conclusion that the X-ray emission observed in 
the 6-10 keV bandpass (encompassing three prominent iron 
lines Fe I Ka, Fe xxv Kq, and Fe xxvi Lya at 6.4, 6.7, and 
6.9 keV respectively plus a very-hard continuum), to a large 
extent follows a smooth, symmetric, centrally-concentrated 
spatial distribution consistent with an origin in unresolved 
point sources. There are two regions where additional emis- 
sion components are clearly required. In the Long +ve re- 
gion. X-ray fluorescence and reflection from dense molecular 
clouds is a major contributor to the 6.4-keV line emission 
and the very-hard continuum band. In this same region, the 
excess 6.7-keV line flux may be interpreted in terms of soft- 
thermal emission, with the lack of an associated excess at 
6.9-keV constraining the temperature of the plasma to be 
less than a few keV. Excess 6.7 keV emission is also seen in 
the Lat -ve region, in a structure which may traced back to 
the central 10 pc region. 

The next step was to investigate whether the spectral 
characteristics of the centrally-concentrated component are 
consistent with a putative origin in the integrated emission 
of low- luminosity sources. An indication of the spectral char- 
acteristics of this component is provided by the source model 
normalisations determined for each of the bands considered 
in Fig. |4l The normalisation of the source model for the 
6.7-keV band relative to that measured in the 7.2-10 keV 
band implies an equivalent width for the 6.7-keV iron line 
of « 630 eV (assuming the underlying continuum has a 7.5- 
keV thermal bremsstrahlung spectral form). Similarly, the 
source model normalisations measured for the 6.4-keV and 
6.9-keV bands imply equivalent widths in excess of 200 eV 
for both lines. This characterisation of the source emission in 
terms of a hard continuum coupled with a trio of prominent 
iron lines is substantiated by the results of spectral fitting 
presented in the next section. 



3.4 Modelling the X-ray spectra 

The X-ray spectra measured by the pn camera in the four 
regions identified in Fig. [2] (which we again refer to as the 
Long +ve. Long -ve, Lat -|-ve, Lat -ve regions) are shown in 



Fig. [6l In the spectral fitting of these data we have consid- 
ered the three emission components discussed earlier ( ^3.ip . 

The first component represents the integrated emission 
from unresolved sources and was characterised in the spec- 
tral fitting as a bremss continuum with kT — 7.5 keV, plus 
Gaussian-line components corresponding to the Fe I Ka, Fe 
xxv Ka, and Fe xxvi Lya lines. The intrinsic width of 
the 6.4-keV and 6.9-keV lines was fixed at 30 eV, in line 
wit h previous estimates of apparent line-widths in pn spec- 
tra (|Capelh et al.ll2012l '). In the case of the 6.7-keV He-like 
iron line, preliminary fitting of the Lat -|-ve spectrum in- 
dicated a best-fit value for the intrinsic width of ~ 55 eV, 
and in the subsequent spectral analysis the intrinsic width of 
this line was fixed at this value. Similarly the energies of the 
three iron lines were measured to be 6.41 ± 0.01 keV, 6.68 
± 0.01 keV, and 6.97 ± 0.01 keV, and subsequently fixed at 
these values. An Fe I K/3 line was included at an energy of 
7.07 keV and width 30 eV with a n ormalisation tied to 11 
per cent of that of the Fe I Ka line l|Kovama et al.ll2009l ). 

The second emission component representing the fiu- 
orescence and reflection from dense molecular clouds was 
modelled as a power-law continuum with photon index F 
plus two Gaussian lines corresponding to neutral Fe Ka and 
K/3 (with parameters other than the normalisation of the Fe 
Ka line set as above). 

Soft thermal emission constituted the final ingredient of 
the composite spectral model. Twin apec components with 
temperatures fixed at fcTi = 0.8 keV and kT2 = 1.5 keV were 
used to represent the likely complex temperature structure. 
In the spectral fitting the metal abundance was allowed to 
vary from region to region, but in all cases was tied across 
the two thermal components. 

It was assumed that each of above components was sub- 
ject to absorption. For the soft thermal emission we set 
A^H ~ 6.5 X lO^'^ cm"'^, in line with recent estimates of 
the s oft X-ray absorption along the line-of-sight to the GC 
(e.g.. IMuuo et al.. 2006. '). In the case of the integrated emis- 
sion of the unresolved sources and the refiection/scattering 
from molecular clouds, the absorption column density was 
increased to Nu = 12 x 10^^ cm~^, so as to include an ad- 
ditional contribution intrinsic either to the source spectra 
or the molecular cloud environment. Unfortunately the in- 
trinsic absorption is rather poorly constrained due to the 
pervasive soft thermal emission, thus a value of 5.5 x 10^^ 
cm"'^ was eventually chosen so as to give reasonably consis- 
tent results across the four regions. 

With all the elements of the composite spectral model in 
plac43, we first attempted to fit the spectra of the Lat +ve 
and Long -ve regions with the integrated source emission 
model (involving four free parameters, i.e., the normalisa- 
tions of the bremsstrahlung continuum and the three iron 
emission lines) together with a kT = 0.8 keV apec com- 
ponent (involving two further free parameters - the apec 
normalisation and metal abundance). This restricted model 
provided a good match to the observational data (Fig. [6j 

^ Our spectra model also included a power-law continuum with 
photon index F = 1.4, subject to an absorbing column density 
A'^H = 13x lO'^'^ cm~'^ w i th a normalisation matche d to that of the 
CXB l lLumb et al.ll2002l : lDe Luca fc Molendill2004l ') . However, the 
inclusion of this component had negligible impact on the spectral 
fitting results. 
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Figure 6. X-ray spectra measured by the pn camera in observation 0202670801 in four regions - Upper-left panel: the Lat -|-ve region; 
Upper-right panel: the Long -ve region; Lower-left panel: the Lat -ve region; Lower-right panel: the Long +ve region. Each panel shows 
the X-ray spectrum after background subtraction in the 2—10 keV band, excluding the region around the Cu Ka instrumental line. 
The best-fitting model and the fitting residuals are also shown. The contribution of different spectral components is illustrated by the 
different colour curves - blue: the 7.5 keV thermal bremsstrahlung continuum plus iron lines associated with unresolved sources; green: 
the fluorescent iron line and reflected continuum from dense molecular clouds; red: soft thermal emission with kT = 0.8 keV; orange: 
(relatively) soft thermal emission with kT = 1.5 keV. 



with best-fitting parameters as listed in Table [S] Of par- 
ticular note is the fact that the derived equivalent widths 
for each line are reasonably consistent across the two re- 
gions, and broadly in line with the expectations based on 
the spatial cuts analysis ( i)3.3p . Based on our best-fitting 
unresolved source model, a 2-10 keV flux of 1.1 x 10~^^ erg 
s~^ cm~^ has, on average, an associated 6.7-keV photon flux 
of 4.0 X lO"'^ photon s"^ cm"^. 

In extending the spectral analysis to the Lat -ve region, 
it was necessary to introduce the hotter kT = 1.5 keV apec 
component into the fltting (thus accounting for the excess 
6.7-keV iron line emission observed in this region - see i]3.3|) . 
The best-fit model for this region is again illustrated in Fig. 
[6] and the best-fit parameter values tabulated in Table (5] 
Again the parameters for the integrated source emission are 
strikingly similar to those of the Lat -|-ve region. 

Finally, for the Long -|-ve region we required the full 
composite spectral model, consistent with our earlier identi- 
fication of this region as the most active of the four consid- 
ered. Here the spectral-fitting approach was to fix the pa- 
rameters of the unresolved-source emission model at those 



obtained from the Long -ve region. The additional free- 
parameters were those quantifying the reflection from the 
molecular clouds (i.e., the continuum normalisation and 
photon index) and the fluorescence (i.e., the normalisation 
of the 6.4-keV line) Again, Fig.|6]and Table[5]summarise the 
results. The derived photon index of F « 1.8 is consistent 
with other recent estimates (|CapeUi et al.ll2012l ). Similarly 
the inferred large equivalent width of the fluoresce nt line 
(~ 1 . 7 keV) is in line with p ublished values ( Ponti et al.l 
I2OIOI : ICapeUi et aLll201ll . |2012| ) . 

An interesting result from Table [5] is the inferred vari- 
ation in the metal abundance across the four regions, with 
values ranging from 0.5-2 Zq. We measure the highest metal 
abundance to the north-east of Sgr A* (Long -l-ve), the 
site with the most recent supernova activity, which could 
have led to a localised enrichment of the thermal plasma. A 
Suzaku study of the GC X-ray emission in the same region 
determined the abund ance to be in the range ~ 1-2 Zq 
(|Nobukawa et all |201G| ) , consistent with the current mea- 
surements. 

In order to test the hypothesis that the integrated 
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Table 5. The best-fitting spectral parameters for the four regions. 



Region 


Lat +ve 


Long -ve 


Lat -ve 


Long +vc 


Unresolved source component 


Nil (xlO^^ cm-2) 
kT (keV) 
Norm. (xlO~^) 
EWeA (eV) 
EWe.7 (eV) 
EW6.9 (eV) 


12 
7.5 
5.90 ±0.09 
236 ± 27 
691 ± 37 
296 ± 33 


12 
7.5 
6.93 ±0.09 
233 ± 23 
781 ± 37 
283 ± 27 


12 
7.5 
5.92 ±0.19 
188 ± 30 
710 ± 58 
366 ± 37 


12 
7.5 
6.93 
233 
781 
283 


Rcficction component 


ATh (x1022 cm-2) 

r 

Norm. (xlO~'') 
EWeA (eV) 








12 

1655 ± 80 


Soft thermal plasma component 


Nil (xlO^^ cm~2) 
kTi (koV) 
Norm. (xl0"2) 
kT2 (keV) 
Norm. (xl0~2) 
Z (Zq) 


6.5 
0.8 

9 6+^ ° 

49+0 " 


6.5 
0.8 

73+0 12 
'J- '•^-0.10 


6.5 

0.8 
^ 7+1.0 

1.5 

r, 7+0.3 
^- ' -0.3 

-0.14 


6.5 
0.8 

1.5 

„ O1+0.24 
U.Sl_o 21 
9 -,0+0.56 




1473.7/1422 


1404.1/1414 


1560.9/1423 


1569.2/1407 



source emission is largely due to magnetic CVs, we have 
employed the spectral model for intermediate polars (IPs) 
developed bv lYuasa et al.l (|2010l ). This IP model has three 
free parameters: the white dwarf (WD) mass Mwd; the 
relative iron abundance Zyc, and a normalisation factor. 
When this spectral model was substituted for the thermal 
bremsstrahlung continuum and the associated lines at 6.7 
keV and 6.9 keV (with the normalisation of the 6.4-keV line 
retained as a free parameter since the IP model does not in- 
clude 6.4-keV iron line emission), we obtained an excellent 
fit to the data with a. /v value for the joint-fitting of the 
spectra for the Lat +ve and Long -ve regions of 2944.8/2849. 
The best-fitting values for the white dwarf mass and iron 
abundance were 0.49 ± 0.02 Mq and 0.97 ± 0.05 Zq. This 
estimate of the typical mass of a magnetic white dwarf in the 
GC is , in fact, very similar to that obtained bv lYuasa et al.l 
l|2012l) in their analysis of Suzaku XIS and HXD/PIN spec- 
tra pertaining to a set of Galactic bulge fields. However, as 
noted by these authors this mass estimate is very likely bi- 
ased downwards by other types of CV, such as dwarf novae 
in quiescence, which have lower representative temperatures. 
Also the lack of spectral information above 10 keV in the cur- 
rent analysis, severely compromises the measurement of the 
effective temperature of the emission (kT > 15 keV), which 
directly impinges on the mass estimate. Nevertheless, the 
analysis does demonstrate that it is very plausible that the 
bulk of the emission seen above 5 keV in the XMM-Newton 
spectra is attributable to an unresolved source population 
dominated by magnetic CVs. 



4 DISCUSSION 

4.1 The low- luminosity X-ray source population 
at the Galactic Centre 

Our spatial and spectral analysis confirms the view that 
much of the hard X-ray emission observed within 20' of Sgr 
A* can be attributed to the integrated emission of unre- 
solved low-luminosity X-ray sources. The spectral character- 
istics of this component - a hard continuum (with kT 7.5 
keV when modelled as thermal bremsstrahlung) plus three 
iron lines (Fe I Ka, Fe xxv Ka, and Fe xxvi Lya) of rel- 
atively high equivalent width - places a tight constraint 
on the nature of the underlying population. These spec- 
tral properties rule-out a significant contribution from most 
classes of coronall y-active stars a nd binaries or from young 
stellar objects (jMuno et al.ll200'3 ). High-mass and low-mass 
X-ray binaries (HMXBs and LMXBs) represent the most 
X-ray lumi nous X-ray p opulation within the Galactic disc 
and bulge (|Muno et al.| [2003: and references therein), but it 
seems highly unlikely that very low-luminosity versions of 
such sources exist in sufficient numbers in the GC to ex- 
plain the observations, although some quies cent LMXBs do 
seem to have the requisite hard spectra fe.g.. lWiinands et al.l 
f2005,). 

Currently, the leading candidate for the unresolved 
source population is magnetic CVs (|Muno et al ] l2004al : 
iLavcock et"alll2005l ). which in principle include both polars 
and IPs. Although the spectra of these subclasses of CV can 
both be characterised as multi-temperature thermal plas- 
mas, the plasma temperatures of polars are lower d ue to en- 
hanced cyclotron cooling (e.o.. lCropper et al1ll998l ) and it is 
the IPs which likel y dominate the har d (> 10 keV) emission 
seen from the GC (|Yuasa et al.ll2012l ) and, correspondingly. 
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make the largest contribution to the key iron-line tracers. 
Nearby, well-studied magnetic cataclysmic variables typi- 
cally have 2-10 keV X-ray luminosities in the range from 
Iq3i _ xo34 erg s"^ (e.g.. iSazonov et aljliooll , although 
the true low-luminosity bound for the population may be 
at least one, and possib ly two, orders of magnitude fainter 
(|Revnivtsev et al]|2009l V A key question for the unresolved 
source model is whether there are sufficient magnetic CVs in 
the GC to account for the bulk of the observed hard X-ray 
luminosity. 

As a starting point, we have determined the integrated 
X-ray luminosity, Lx, of the unresolved source component. 
For this purpose, we first measured the normalisation of the 
empirical surface brightness model, when fitted to the Long 
-ve and Lat -f ve segments of the spatial cuts pertaining to 
the 4.5-6 keV pn images (noting that in these regions, the 
unresolved sources appear to provide the dominant contri- 
bution to the 4.5-6 keV flux). Next we converted pn count 
s"^ in the 4.5-6 keV band to X-ray flux in the full 2-10 
keV band (via the spectral model of the unresolved source 
component discussed previously), and applied the necessary 
scaling factors to transform the surface brightness model 
units to X-ray luminosity per square parsec. Finally, we de- 
termined the radial surface-brightness distribution from an 
image derived from the modeO The result is shown in the 
upper-left panel of Fig[71 Integrating this distribution from 
2'-20' (which corresponds to a projected offset of 4.6-46 pc 
at the GC), we obtain Lx ~ 2.0 x 10^® erg s^^ as the total 
integrated (2-10 keV) X-ray luminosity of the unresolved 
sources (below a source detection threshold of Lx ~ 10^'^ 
erg s"^). 

The next step was to consider the GC stellar mass 
contained within this region. Here we u se the same mass- 
model prescription as iMuno et al.l (|2006l ). who in turn use 
parameters obta i ned f rom an analysis of NIR maps by 
iLaunhardt et al.l (|2002l fl 

Within the central few hundred parsecs of the Galaxy, 
two major stellar distributions have been identified: the Nu- 
clear Stellar Cluster (NSC) and the Nuclear Stellar Disc 
(NSD). 

The NSC has a mass density, which varies with radius 
R as: 

^ " 1 + (7?/i?,c)"i ■ 

For i? < 6 pc, Pc = 3.3 x 10*^ M0pc"^ Rc = 0.22 pc, and 
ni = 2. For 6 < ii < 200 pc, ni = 3, Rc remains the same 
and pc must be adjusted so that the profile is continuous at 
-R = 6 pc. 

The mass distribution of the NSD can be modelled as 
a power-law plus exponential function. 



* We use circular annuli of width 1'. Hence the derived radial 
d istribution re p resent s an average over latitude and longitude. 
^ iMuno et alj ll2006t) adopt the same parameters values as 
ILaunhardt et al. I ll2002h despite the fact that the original anal- 
ysi s was based on a GC d istance of 8.5 kpc. However, as noted 
by ILaunhardt et aL I | |2002|'| . in practice, such scaling corrections 
are typically much smaller than the uncertainties in the param- 
eter determinations. In the present pap e r we employ the same 
mass- model formulation as iMuno et alj | |2006|) , along with the 
assumption of a GC distance of 8 kpc. 



p = PdR exp{-\z\/za,). (3) 

For R < 120 pc, pd = 300 Mqpc'^, n2 = 0.1, and Zd = 45 
pc. For 120 < R < 220 pc, becomes n2 = 3.5 with Zd un- 
changed; Pd must be adjusted so that the function is con- 
tinuous for n2 = 3.5. Beyond 220 pc this disc component is 
assumed to cut-off rapidly. 

The other com ponents of the mass-model discussed by 
iMuno et al.l l|2006l ) are the Galactic bulge and disc. However, 
these are not considered here, since they contribute only to 
the constant term in our empirical surface brightness model 
(on the angular scales considered), which was not included 
in the derivation of the integrated luminosity. 

We integrated the mass distributions in the NSC and 
NSD components along a set of lines-of-sight to give the stel- 
lar mass per square parsec on a 2-d grid of positions centred 
on Sgr A*. From this 2-d distribution we then determined 
the average stellar-mass density within circular annuli of 
width 1' as a function of radial offset in the range 2'-20'. 
The result is shown in the upper-right panel of Fig[7l which 
displays both the individual contributions of the NSC and 
NSD and their sum. It is evident that the NSD provides the 
dominant mass contribution for offset angles greater than a 
few arcmin and that the NSC component becomes insignif- 
icant outside ~10'. 

The scale h e ight for the NSD of 45 pc reported by 
iLaunhardt et al.l (|2002l l. based on IRAS 100 [im dust mea- 
surements, corresponds to 18.2' at the distance of the GC 
assumed by these authors. This is in excellent agreement 
with the scale height inferred from the X-ray observations 
(see Fig. [5]) and consistent with the fact that the NSD is 
the major contributor to the stellar mass within the region 
sampled by the X-ray observations. 

The X-ray emissivity per unit stellar mass can be read- 
ily determined from the ratio of the curves presented in the 
upper-left and upper-right (total mass) panels of Fig. [7] The 
result is shown in the lower-left panel of Fig. [7] A smooth 
trend is apparent with offset angle, with the inferred emis- 
sivity varying from « 5 x 10^^ erg s~^ Mq~^ at 20' up to 
almost twice this value at 2\ 

By way of comparison lSazonov et al.l (|2006l ) constructed 
an X-ray luminosity function for faint sources in the local 
neighbourhood, from which they determined the total X- 
ray emissivity for CVs and active binaries (that is sources 
associated with the old stellar population) to be ~ (3.1 ± 
0.8) X 10^'' erg s"^ M©"^ in the 2-10 keV band, with CVs 
contributing roughly one-third of this total. More recently, 
using a deep RXTE / PCA scan across the G alactic Plane 
at III ~ +18.0° . iRevnivtsev fc Molkovl (|2012l ) calculate the 
2-10 keV emissivity in the GRXE to be (3.0±0.3) x 10^'' erg 
s~^ Mq~^. Moving closer to the GC, from a Chandra deep- 
field observation targeted at the inner Ga l actic bulge {In 
- 0.113° , 6ii - -1.424° ). iRevnivtsev et all (|201lh estimate 
the 2-10 keV X-ray emissivity to be (2.4 ± 0.4) x 10^'^ erg 
s~^ M0~'^. Finally from a similar Chand ra study based on 
deep-fi eld data taken within l'-4'of Sgr A*. lRevnivtsev et al.l 
(20071) derived an X-ray emissivity value of ~ (8.5 ± 4.3) x 
10^^ erg s"'^ M©"^ (corrected to a GC distance of 8 kpc). 
Here the large error stems from uncertainty in the GC mass 
model - see the discussion below. 

Taken as a whole, the published measurements provide 
evidence of an up- lift in the X-ray emissivity per unit stellar 
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Figure 7. Upper-left panel: Variation of the 2—10 keV X-ray surface brightness (in units of 10'^^ erg s~^ pc~^) as a function of the 
angular offset from Sgr A*. The thick and thin black lines are the luminosity per square parsec derived from the X-ray images and 
predicted from the mass-model respectively. Upper-right panel: Projected mass density (in units of 10* Mq pc~^) in the GC showing the 
contribution of the NSC (red), the NSD (blue), and their sum (black). Lower-left panel: The inferred X-ray emissivity per unit stellar 
mass (in units of 10^^ erg s~^ Mq~^) as a function of the angular offset from Sgr A*. The dotted line is the "base-level" X-ray emissivity 
which is representative of measured GRXE values. Lower-right panel: The variation in the stellar mass density on the plane (in units of 
M0 pc~^) for the GC region as a function of the d istance (in pc) from Sgr A*. The solid black curve represents the sum of the NSC 
(red) and the NSD (blue) components in the lMurio ^ ct al. (2006) model. The dotted black curve is the mass-model (sum of the NSC and 
NSD components) inferred on the basis of the assumption that the X-ray emissivity in the region is 3 X lO^'^ erg s~^ M©"^. 



mass of the unresolved source population within < 1 deg of 
Sgr A*, which is substantiated by our current results. 

Two interpretations are possible for the upward trend 
in the X-ray emissivity per unit stellar mass within the cen- 
tral 100-pc zone evident in Fig. [T] The first is that some 
characteristic of the underlying source population is chang- 
ing in tandem with the stellar density. If we assume the 
"base-level" X-ray emissivity of the old stellar population 
is roughly 3.0 x lO^'^ erg s~^ Mq~^ (i.e., a value repre- 
sentative of the GRXE excluding the GC) and apply this 
scaling to the projected mass distribution (Fig. [T] top-left 
panel), then the observed X-ray luminosity shows a factor 
~ 2 enhancement relative to this base-level. Conceivably 
the binarity of stellar population, and hence the incidence 
of magnetic CVs, might be dependent on the stellar density. 
Alternatively, the contribution from other source types {e.g., 
extreme RSCVn binaries) might rise within 50 pc of Sgr A* 
in a f ashion analogous to t hat inferred for the central stellar 
cusp (jSazonov et al.ll2012l ). The nucleus of M31 presents an 



interesting comparator. Within 1' of the centre of M31, the 
incidence of luminous X-ray sources (Lx < 10'^^ erg s~^) 
increases more rapidly than might be predicted from the K- 
band Ught l|Voss fc Gilfanovll2007ah . a result which has been 
interpreted in terms of the dynami cal formation of LMXB s 
within a dense stellar environment (|Voss fc Gilfanovll2007bl ) . 
At the distance of M31, 1' corresponds to a linear extent of 
200 pc, a scale-size not dissimilar to the extent of the NSC 
and NSD components discussed above. Plausibly, dynamical 
processes may have given rise to a higher incidence of X-ray 
emitting CVs and active binaries within the central 100-pc 
region, than that pertaining in the Galaxy as a whole. 

An alternative explanation of the observed upward 
trend in the X-ray emissivity per unit stellar mass and the 
implied enhanced X-ray luminosity of the central 100-pc 
region, is that the mass model we have employed is inac- 
curate. In fact a number of authors have noted that the 
uncertainties in the current mass model, including the rel- 
ative scaling of the NSC and NSD components, may be as 
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high as a factor of ~ 2 llLaunhardt et al.ll2002l : iMuno et al.l 
I2OO6I : iRevnivtsev et al.ll2007l '). 

We have investigated what changes in the mass distri- 
bution are imphed by our measurements, if we make the 
assumption that the X-ray emissivity within the GC re- 
gion is constant. Bearing in mind that we are deahng with 
quantities measured in projection on the plane of the sky 
(i.e., integrated along the line-of-sight through the whole 
nuclear region), we have necessarily sought a representative 
solution rather than a unique one. In that context, we find 
that the upward trend in the emissivity can be counteracted 
by steepening, the power-law function describing the de- 
crease with radius of the NSD mass density within R < 120 
pc; specifically the power-law index must be changed from 
n2 = 0.1 to 712 ~ 1.2. The total GC mass profile which re- 
sults in a roughly constant X-ray emissivity of 3.0 x 10^^ 
erg s~^ M0~^ between 2'-20', is shown in the lower-right 
panel o f Fig. \7\ We note that compared to the c urrent mass 
model ijLaunhardt et al.ll2002l : lMuno et all2006l 'l. additional 
mass (roughly a factor 2) is required within 10-100 pc of Sgr 
A*; however, it is unclear whether this is appropriately as- 
signed to the NSG, the NSD or to a combination of both. 
In summary, this analysis demonstrates that X-ray tracers 
may eventually provide a very effective way of constraining 
the stellar mass distribution in the GC; of course, this is 
subject to the caveat that we need to better understand the 
nature of GC X-ray source population and the evolutionary 
influences to which it has been subject. 

Finally we may use the estimates of this section to illus- 
trate just how extreme the environment is at the centre of 
the Galaxy. If we take the stellar mass density to be ~ 10^ 
Mq pc~"^ (the value at ~ 20 pc on the plane in our revised 
mass model), then the base- level X-ray emissivity per unit 
stellar mass assumed above implies an X-ray volume emissiv- 
ity of 3 X 10'^" erg s~^ pc~^. If a typical magnetic CV in this 
region has an X-ray luminosity of Lx ~ 10'^" erg s~^, con- 
sistent with the faint end of the so urce distribution probe d 
in Chandra deep-field observations (|Revnivtsev et al1l2009f ). 
then the volume density of magnetic CVs is ~ 3 pc~^. How- 
ever, only ~ 5 per cent of all CVs are thought to be IPs 
ijRuiter et al.|[200^ ). so the projected number density of all 
types of CVs rises to 60 pc~^, compared to the inferred local 
densi ty of just ~ lO""^ pc"^ l|Warneilll995l : ISchwope eFsI} 
|2002| ). 

4.2 Absence of very-hot diffuse thermal plasma in 
the GC 

The GRXE can be traced along the Galactic Plane 
from I in I ~ 90° through to t he GC, where its surface 
brightness reaches a m aximum ama et al ] |l986l . ll989l : 
IRevnivtsev et all l2006t ). In recent years the debate as to 
whether this feature is due predominantly to the presence 
of a truly diffuse very-hot (kT ~ 7.5 keV) thermal plasma 
or due to the integrated emission of point sources, has inex- 
orably swung in favour of the latter interpretation. 

A key recent result in the context of the origin of the 
emission from the GC was the discovery that in a deep Chan- 
dra observation at {lu,bu = +0.08° , —1.42°), more than 
~ 80 per cent of the continuum and line e mission around 6.7 
keV could be resolved into point sources (|Revnivtsev et al] 
I2OO9I ). Also detailed spectral analysis of Suzaku observations 



in the Galactic bulge region, suggests that virtually all of the 
unresolved X-ray emission can be ex plained in terms o f the 
integrated emission of point sources (|Yuasa et al.ll2012l l. 

In this paper, we have investigated the distribution of X- 
ray emission in the region between 2' and 20' from Sgr A* in 
a variety of bands, encompassing three iron-line tracers and 
also the very-hard 7.2-10 keV continuum. We find that in 
many regions the emission in these bands follows a smooth, 
symmetric centrally-concentrated spatial distribution. Fur- 
thermore, to within a factor of 2, the observed surface bright- 
ness distribution matches the projected mass distribution in 
the old stellar population in the G C, speciflcally that in 
the Nuclear Stella r Disc component (|Launhardt et al"]|2002l : 
iMuno et al.ll200^ ). We interpret this as strong support for 
the conjecture that the observed emission is predominantly 
the integrated emission of unresolved point sources. 

Excluding the central 2', where the Sgr A East SNR, 
the Central Cluster, and other features combine to produce 
a com plex morphology (|Baganoff et al. I I2003I : M 

uno et al.l 

l2004bl ). and also a few regions dominated by individual X- 
ray binaries or other bright sources (e.g., the Arches Clus- 
ter), we identify two regions within 20' of the GC where 
the underlying glow of point sources is confused by ad- 
ditional emission components. The first is the north-east 
region where X-ray fiuorescence and refiection from dense 
molecular clouds gives rise to a strong signal in the 6.4-keV 
iron-line and the very-hard continuum bands. Excess emis- 
sion at 6.7-keV can be seen in this north-east region and, 
also in a second "anomalous" region to the south of Sgr A*. 
In both cases we interpret the 6.7-keV excess as a concentra- 
tion of soft thermal emission with temperature ~ 1.5 keV. 
The key point is that within this overall setting there is no 
requirement for any excess very-hard thermal emission which 
might be attributed to a residual diffuse GRXE component. 

In a recent paper, lUchivama et al.l |20l3) used Suzaku 
measurements to trace the distribution of iron-K emission 
from the GC proper through to regions of the inner Galac- 
tic Plane {\lu\ = 1° — 10° ). Similar to the current anal- 
ysis, these authors compared the spatial variation of the 
Fe XXV Ka line with a stellar mass distribution model. 
They found that when the X-ray emissivity per unit stel- 
lar mass was fixed at the value required to explain the 
observed emission outside of \lii\ ~ 1.5° , then the X-ray 
emission observed within 1° was 3.8 times under-predicted. 
They concluded that either a new source population with 
extremely strong Fe xxv Ka emission was required in the 
GC or that the majority of hard GC X-ray emission must 
originate in a truly diffuse optically-thin thermal plasma. 
In the present paper we argue that the X-ray emissivity in 
the GC is roughly a factor 2 higher than esti mates pertain- 
ing t o the local Galajcy and Gala ctic Plane (jSazonov et al.l 
I2OO6I : IRevnivtsev fc Molkovl [201^ ') . A possible explanation 
of these different estimates of the degree to which the GC 
X-ray emissivity is enh anced (3.8x as opposed to 2x) is that 
lUchiv ama et all l|201ll ) use a scaling for their Galactic disc 
component, taken from iMuno et all (|2006l ). which is based 
on a total mass in the G alactic disc of 10^^ Mq. However, 
as noted bv iHong et al.l (2009) this is likely to be an over- 
estimate of t he to tal stellar mass in the disc; for example 
iRobin et al.] l|2003f) quote a value of 2.2 x 10^° Mq. What- 
ever the actual degree of enhancement of the X-ray emissiv- 
ity in the central 100-pc region, we hold to the view that 
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the observed spatial distribution and spectral properties of 
this component are best explained in terms of the integrated 
emission of unresolved sources. 



5 CONCLUSIONS 

We have used the extensive set of observations available in 
the XMM-Newton archive to make mosaiced images of the 
central 100 pc x 100 pc region of our Galaxy in a variety of 
broad and narrow bandpasses within the 2-10 keV energy 
range. 

Using three iron-K lines (Fe I Ka at 6.4 keV, Fe xxv Ka 
at 6.7 keV and Fe xxvi Lya at 6.9 keV) as tracers, we decom- 
pose the GC X-ray emission into three distinct components. 
The first of these components is the integrated emission of 
unresolved point sources with 2-10 keV X-ray luminosity 
Lx < 10"^"^ erg s~^. The second is the fluorescent line emis- 
sion and the reflected continuum from dense molecular ma- 
terial. The third is the soft diffuse emission from thermal 
plasma in the temperature range, kT « 0.8-1.5 keV. We 
show that the spatial distribution and spectral properties of 
the GG X-ray emission (excluding the central 2' region and 
a few other regions confused by luminous sources) can, to 
a good approximation, be described in terms of these three 
components. 

A key finding is that unresolved sources account for 
much of the observed emission in 6-10 keV band, including 
the iron-line emission. We find that the surface brightness 
of this component falls off as ~ 6~^'^ , where 9 is the angular 
offset from Sgr A*, coupled with a latitudinal scale height 
factor (^sc=18.6ti'2 arcmin. 

We assume that the unresolved sources are associ- 
ated with the old stellar population of the Galaxy, which 
can be traced in the GG through the NIR light, on 
th e basis of which, a 3-d mass model has been derived 
bv iLaunhardt et al.1 (|2002l ). According to ILaunhardt et al] 
(|2002l ) there are two major stellar structures within the cen- 
tral 100 pc of the Galaxy, namely the Nuclear Stellar Cluster 
(NSG) and the Nuclear Stellar Disc (NSD). When we project 
these two components of the mass-model onto the 2-d plane 
of sky, we find that it is the NSD which dominates the re- 
gion sampled by the X-ray measurements. The X-ray surface 
brightness distribution (as parameterised above) rises more 
rapidly with decreasing 6 than the mass-model prediction. 
One interpretation is that the 2-10 keV X-ray emissivity 
increases from « 5 x 10^^ erg s~^M0~^ at 20' up to al- 
most twice this value at 2'. Alternatively, some refinement 
of the mass model may be required. Specifically if we as- 
sume the X-ray emissivity in the central region is similar to 
that pertaining to GRXE as a whole (i.e., ~ 3 x 10^^ erg 
s~^M0~^), then roughly a factor of two additional mass is 
required within 10-100 pc of Sgr A*. However, it is unclear 
whether this applies to the NSG, the NSD, or a combination 
of both components. 

When we set 712 = 1-2 in the NSD model and revise 
the model normalisation so as to conserve the total mass 
in the NSD within our 2'-20' field of view, we obtain an 



X-ray emissitivity per unit stellar mass of « 6.5 x 10 erg 
s~^M0~^. This is within a factor 2 of other estimates of this 
quantity based on the local X-ray luminosity function and 
studies of the GRXE component outside of the GG region. 



The unresolved hard X-ray emitting source population, 
on the basis of spectral comparisons, is most likely dom- 
inated by magnetic cataclysmic variables, primarily inter- 
mediate polars. By way of example, we use the X-ray ob- 
servations to estimate that there are 60 GVs (of all types) 
per cubic parsec at a radial distance of 20 pc from Sgr A*. 
This is close to seven orders of magnitude over the local GV 
number density. 

A further major finding of our work is that our com- 
posite emission model, encompassing the three components 
described above, provides a sufficiently good description of 
the observations so as to rule out any substantial additional 
source of hard X-ray luminosity in the GG. It has long been 
conjectured that a significant fraction of the hard X-ray 
emission from the GG originates in very-hot (~ 7.5 keV) dif- 
fuse thermal plasma. Within our modelling framework, the 
only way of accommodating a substantial contribution of 
this nature, is within the emission component we attribute 
to unresolved sources. However, the properties of the lat- 
ter, namely the smooth, symmetrical, centrally-concentrated 
spatial distribution, an X-ray emissivity reasonably in ac- 
cord with local estimates, and an X-ray spectrum matching 
that of intermediate polars, would seem to exclude the dif- 
fuse, very-hot thermal plasma hypothesis as a viable option. 
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